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Abstract Titan’s polar surface is dotted with hundreds of lacustrine depressions. 
Based on the hypothesis that they are karstic in origin, we aim at determining the efficiency 
of surface dissolution as a landshaping process on Titan, in a comparative planetology 
perspective with the Earth as reference. Our approach is based on the calculation of 
solutional denudation rates and allow inference of formation timescales for topographic 
depressions developed by chemical erosion on both planetary bodies. The model depends 
on the solubility of solids in liquids, the density of solids and liquids, and the average annual 
net rainfall rates. We compute and compare the denudation rates of pure solid organics 
in liquid hydrocarbons and of minerals in liquid water over Titan and Earth timescales. 

We then investigate the denudation rates of a superficial organic layer in liquid methane 
over one Titan year. At this timescale, such a layer on Titan would behave like salts or 
carbonates on Earth depending on its composition, which means that dissolution processes 
would likely occur, but would be 30 times slower on Titan compared to the Earth due to 
the seasonality of precipitation. Assuming an average depth of 100 m for Titan’s lacustrine 
depressions, these could have developed in a few tens of millions of years at polar latitudes 
higher than 70° N and S, and a few hundreds of million years at lower polar latitudes. The 
ages determined are consistent with the youth of the surface (< 1 Gyr) and the repartition 
of dissolution-related landforms on Titan. 
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1 Introduction 


Along with the Earth, Saturn’s icy moon Titan is the only planetary body of the entire 
Solar System that possesses lakes and seas {Lopes et ai, 2007; Stofan et ai, 2007; Hayes 
et ai, 2008). Some of these lakes and seas are currently covered by liquids, while others 
are not {Hayes et al, 2008). Most are located in the polar regions {Hayes et al., 2008; 
Aharonson et al, 2009), although a few occurrences have been reported at lower latitudes 
{Moore and Howard, 2010; Vixie et ai, 2012). The currently-filled lakes and seas are 
located poleward of 70° of latitude in both hemispheres whereas most empty depressions 
are located at lower latitudes (Figure 1) {Hayes et ai, 2008; Aharonson et al., 2009). 

Altogether, empty depressions, lakes, seas and fluvial channels, argue for the presence 
of an active “hydrological” cycle on Titan similar to that of the Earth, with exchanges 
between the subsurface (ground liquids), the surface (lakes, seas, fluvial channels) and Ti¬ 
tan’s methane-rich atmosphere, where convective clouds and sporadic intense rainstorms 
have been imaged by the Cassini spacecraft instruments {Turtle et al., 2011a). Methane, 
rather than water as on Earth, probably dominates the cycle on Titan {Lunine et al., 2008) 
and thus constitutes one of the main components of the surface liquid bodies observed in 
the polar regions {Glein and Shock, 2013; Tan et al., 2013). Ethane, the main photodis¬ 
sociation product of methane {Atreya, 2007), is also implied in Titan’s lakes chemistry, as 
predicted by several thermodynamical models {Lunine et al, 1983; Raulin, 1987; Dubouloz 
et al., 1989; Cordier et al., 2009; Tan et al, 2013; Cordier et al., 2013a; Glein and Shock, 
2013), or as identified in Ontario Lacus thanks to the Cassini/VIMS instrument {Brown 
et ai, 2008). 

Titan’s lakes are located in topographic depressions carved into the ground by geologi¬ 
cal processes that are poorly understood to date. The origin of the liquid would be related 
to precipitation, surface runoff and underground circulation, leading to the accumulation 
of liquids in local topographic depressions. 

In the present work, we aim to constrain the origin and the age of these depressions. 
Section 2 first provides a brief overview of their geology and a discussion of their possible 
origin based on their morphological characteristics and from considerations about Titan’s 
surface composition and climate. Based on this discussion, we propose a new quantitative 
model, whereby the depressions have formed by the dissolution of a surface geological layer 
over geological timescales, such as in karstic landscapes on Earth. 

In terrestrial karstic landscapes, the maximum quantity of mineral that can be dis¬ 
solved per year, namely the solutional denudation rate, can be computed using a simple 
thermodynamics-climatic model presented in Section 3. The denudation rate depends 
on the nature of the surface material (solubility and density of the minerals) and on 
the climate conditions (precipitation, evaporation, surface temperature). Using this sim¬ 
ple model, it is possible to determine theoretical timescales for the formation of specific 
karstic landforms on Earth, which are compared to relative or absolute age determinations 
in Section 4. 

We apply the same model to Titan’s surface in Section 5. Section 5.1 is dedicated to 
the comparative study of denudation rates of pure solid organics in pure liquid methane, 
ethane and propane, and of common soluble minerals (halite, gypsum, anhydrite, calcite 
and dolomite), cornerstones of karstic landscapes development on Earth, in liqnid water 
over terrestrial timescales. Section 5.2 describes the computation of denudation rates of 
pure solids and of mixed organic surface layers in liquid methane over Titan timescales by 
using the methane precipitation rates extracted from the GCM of Schneider et al. (2012). 
Based on these denndation rates, we compute the timescales needed to develop the typical 
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100 m-deep topographic depressions observed in the polar regions of Titan (Section 6) and 
compare them to timescales estimated from other observations (e.g. crater counting, dune 
formation). 

2 Geology of Titan’s lacustrine depressions 

2.1 Geomorphological settings 

Seas and lacustrine depressions strongly differ in shape (Figure 2). On one hand, seas 
are large (several hundred km in width) and deep (from 150 to 300 - 400 m in depth) 
{Lorenz et al., 2008, 2014; Mastrogiuseppe et al, 2014). They possess dendritic contours 
and are connected to fluvial channels (e.g. Ligeia Mare, Figure 2a) {Stofan et al, 2007; 
Sotin et al., 2012; Wasiak et al, 2013). They seem to develop in areas associated with 
reliefs, which constitute some parts of their coastlines. 

On the other hand. Titan’s lacustrine depressions (Figure 2b) develop in relatively flat 
areas. They lie between 300 and 800 meters above the level of the northern seas {Stiles 
et al, 2009; Kirk et al, 2012). They are typically rounded or lobate in shape and some 
of them seem to be interconnected {Bourgeois et al., 2008). Their widths vary from a few 
tens of km, such as for most of Titan’s lacunae, up to a few hundred km, such as Ontario 
Lacus or Jingpo Lacus. Their depths have been tentatively estimated to range from a 
few meters to 100 - 300 meters {Hayes et al., 2008; Kirk and Howington-Kraus, 2008; 
Stiles et al, 2009; Lorenz et al, 2013), with “steep”-sided walls {Mitchell et al, 2007; 
Bourgeois et al., 2008; Kirk and Howington-Kraus , 2008; Hayes et al., 2008). The liquid- 
covered depressions would lie 250 meters below the floor of the empty depressions {Kirk 
et al., 2007), which could be indicative of the presence of an alkanofer in the sub-surface, 
analog to terrestrial aquifers, filling or not the depressions depending on their base level 
{Hayes et al., 2008; Cornet et al., 2012). The depressions sometimes possess a raised rim, 
ranging from a few hundred meters up to 600 meters in height {Kirk et al., 2007; Kirk and 
Howington-Kraus , 2008). All these numbers are likely subject to modification following 
future improvements in depth-deriving techniques. 

2.2 Geological origin of the depressions 

The geological origin of the topographic depressions and how they are fed by liquids are still 
debated. The geometric analysis of the lakes by Sharma and Byrne (2010, 2011) led to the 
conclusion that, unlike on Earth, the formation mechanism of the lacustrine depressions 
cannot be derived from the analysis of their coastline shapes. Recently, Black et al. (2012) 
and Tewelde et al. (2013) showed that mechanical erosion due to fluvial activity would have 
a minor influence in landscape evolution on Titan. Given this context, several hypotheses 
are being explored to understand how Titan’s lakes have formed. These include: 

1. cryovolcanic origin {Mitchell et al, 2007; Wood et al, 2007), forming topographic 
depressions in which lakes can exist, such as in terrestrial calderas {Acocella, 2007) 
or maars {Lorenz, 1986); 

2. thermokarstic origin {Kargel et al, 2007; Mitchell et al, 2007; Harrisson, 2012), 
where the cyclic destabilization of a methane frozen ground would form topographic 
depressions, such as in periglacial areas on Earth where the permafrost cyclically 
freezes and thaws and forms thermokarst lakes, pingos or alases {French, 2007) ; 
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3. solutional origin [Mitchell et al., 2007; Bourgeois et al., 2008; Mitchell, 2008; Mitchell 
and Malaska, 2011; Malaska et al, 2011; Barnes et al., 2011; Cornet et al, 2012), 
where processes analogous to terrestrial karstic dissolution create topographic de¬ 
pressions, such as terrestrial sinkholes/dolines, playas and pans under various cli¬ 
mates [Shaw and Thomas, 2000; Ford and Williams, 2007). 

On one hand, the general lack of unequivocal cryovolcanic features on Titan tends 
to limit the likelihood of the cryovolcanic hypothesis [Moore and Pappalardo, 2011). A 
methane-based permafrost would be difficult to form on Titan due to the presence of 
nitrogen in the atmosphere [Lorenz and Lunine, 2002; Heintz and Bich, 2009). Its putative 
cyclic destabilization would also be challenging, given the tiny temperature variations 
between summer and winter, day and night, equator and poles [Jennings et al., 2009; Lora 
et al., 2011; Cottini et al., 2012) over all timescales [Aharonson et al, 2009; Lora et al., 
2011). On the other hand, solid organics have been shown to be quite soluble in liquid 
hydrocabons under Titan’s surface conditions [Lunine et al., 1983; Raulin, 1987; Dubouloz 
et al., 1989; Cordier et al, 2009, 2013a,b; Glein and Shock, 2013; Tan et al, 2013), provided 
that they are available at the surface. The observation of bright terrain around present 
lakes and inside of empty depressions, analogues of terrestrial evaporites produced by 
the evaporitic crystallization of dissolved solids, also strengthens this hypothesis [Barnes 
et al, 2011; MacKenzie et al., 2014). 

On Earth, dissolution-related landforms are not restricted to sinkholes/dolines, pans or 
playa, which characterize relatively young karsts [Ford and Williams, 2007). Spectacular 
instances of reliefs nibbled by dissolution, known as cone/cockpit karsts, fluvio karsts 
or tower karsts, exist under temperate to tropical/equatorial climates, such as in China 
[Xuewen and Weihai, 2006; Waltham, 2008), Indonesia [Ford and Williams, 2007) or 
the Carribeans [Fleurant et al., 2008; Lyew-Ayee, 2010). The observation of possible 
mature karst-like terrains in Sikun Labyrinthus (Figure 2c) by Malaska et al. (2010), 
similar to these terrestrial karstic landforms also gives further credence to the hypothesis 
that lacustrine depressions on Titan are karstic in origin. 

2.3 Composition of Titan’s solid surface 

Titan’s surface can be divided into five main spectral units identified by the Cassini/VIMS 
instrument: bright terrain, dark equatorial dune fields or dark-brown units, blue units, 
5 /im-bright units and the dark lakes [Barnes et al, 2007; Stephan et al., 2009). In the 
polar regions, the solid surface appears dominantly as bright terrain, lakes, and patches 
identified as the 5 ^m-bright unit in VIMS data [Barnes et al., 2011; Sotin et al., 2012; 
MacKenzie et al, 2014). The spectral characteristics of the VIMS 5 /xm-bright unit seen 
inside and around some polar (and equatorial) lacustrine depressions indicate the presence 
of various hydrocarbons and nitriles [Clark et al., 2010; Moriconi et al., 2010) and are not 
compatible with the presence of water ice [Barnes et al., 2009). 

The origin of the organic materials is probably linked to the atmospheric photochem¬ 
istry, which results in the formation of various hydrocarbons and nitriles [Lavvas et al., 
2008a,b; Krasnopolsky , 2009) detected by Cassini [Cui et al, 2009; Magee et al, 2009; 
Clark et al., 2010; Coustenis et al., 2010; Vinatier et al., 2010; Cottini et al., 2012). Table 
1 gives the estimated fluxes of some produced organics, as derived from several models. 
Most of these compounds could condense as solids and sediment onto the surface over ge¬ 
ological timescales [Atreya, 2007; Malaska et al., 2011). Most of them would be relatively 
soluble in liquid alkanes [Raulin, 1987; Dubouloz et al., 1989; Cordier et al, 2009, 2013a; 
Clein and Shock, 2013; Tan et al., 2013). 
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It is therefore reasonable to assume that a superficial soluble layer, composed of or¬ 
ganic products, exists at the surface of Titan. Episodic dissolution of this layer would be 
responsible for the development of karst-like depressions and labyrinthic terrains {Bour¬ 
geois et al, 2008; Malaska et al, 2010, 2011; Mitchell and Malaska, 2011; Cornet et ai, 
2012). Evaporitic crystallization could also occur after episodes of dissolution in the liq¬ 
uids, forming evaporite-like deposits {Barnes et al., 2011; Cornet et al., 2012; Cordier 
et al., 2013b; MaeKenzie et al, 2014). 

3 Solutional denudation rates on Earth 

On Earth, karstic landforms develop thanks to the dissolution of carbonate (calcite, 
dolomite) and evaporite (gypsum, anhydrite, halite) minerals under the action of ground- 
water and rainfall percolating through pore space and fractures present in rocks. The 
mineral solubilities vary as a function of the environmental conditions (amount of rain 
and partial pressure of carbon dioxide) {Ford and Williams, 2007). Karstic landforms like 
dolines or sinkholes are often located under temperate to humid climates in carbonates 
(though gypsum or halite karsts also exist on Earth). They reach depths of up to a few 
hundred meters {Ford and Williams, 2007). Karsto-evaporitic landforms like pans are 
located under semi-arid to arid climates. They reach depths of up to a few tens of meters 
{Coudie and Wells, 1995; Bowen and Johnson, 2012). Evaporitic landforms like playas 
are located under arid climates. They are characterized by their extreme flatness and can 
occur in any kind of topographic depression {Shaw and Thomas, 2000). 

Denudation rates (hereafter DR) in terrestrial karstic landscapes are primarily con¬ 
strained by geological (age and physico-chemical nature of rocks) and climate (net pre¬ 
cipitation rates evolution) analyses. Eor many limestone-dominated areas (the majority 
of karst areas), it is commonly assumed that dissolution features are mainly created in 
the epikarstic zone located in the top few meters below the surface {White, 1984; Ford 
and Williams, 2007). Classically, in karstic terrains, the chemical/solutional denudation 
rate DR (in meters per Earth year, or m/Eyr) is related to the rock physico-chemical 
properties and the climate by the following equation {White, 1984; Tucker et al, 2001; 
Ford and Williams, 2007; Fleurant et al, 2008): 

DR = Py, r mca, [m/Eyr] (1) 

Pcalcite 

where is the mass density of liquid water (~ 1000 kg/m^), Mcaicite is the molar 
mass of calcite (in kg/mol), Pcaicite is the mass density of calcite (in kg/m^), r is the mean 
annual net precipitation rate (in m/yr), equivalent to the sum of runoff and infiltration or 
the difference between precipitation and evapotranspiration over long timescales {White, 
2012), and mca is the equilibrium molality of calcite (Ca^"*" cations, in mol/kg), assuming 
an instantaneous dissolution. Eollowing this equation, denudation rates depend linearly 
on the climate precipitation regime and on the molality of dissolved materials. 

Molality calculations are provided in Appendix B for various soluble minerals, based 
on two thermodynamic hypotheses: an Ideal Solution Theory (1ST, no preferential inter¬ 
actions between molecules) and an Electrolyte Solution Theory (EST, preferential inter¬ 
actions between molecules). We incorporated the effect of CO 2 gas dissolved in water, 
which acidifies water due to a series of intermediate reactions (one of which producing 
carbonic acid), and increases the dissolution rates of carbonates (calcite and dolomite, see 
Appendix B for more details). The partial pressure of CO 2 (hereafter Pco 2 ^ m atm) under 
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different climates can be computed as follows {Ford and Williams, 2007): 

logPcOa = -3.47 + 2.09 (1 - ^2) 

where AET is the mean annual evapotranspiration rate in mm/Eyr. Applying the 
above formula gives Pco2 < 3 matm for arid areas, 3 < Pco2 < 12 matm for temperate 
areas and Pco 2 > 12 matm for tropical or equatorial climates. The resulting typical 
values of denudation rates on Earth for common soluble minerals in liquid water are given 
in Table 2. Depending on the climatic conditions, denudation rates of salts are on the 
order of a few hundred /rm to a few cm per Earth year, whereas denudation rates of 
carbonates are on the order of a few /rm to a few hundred /rm per Earth year. Note that 
chemical erosion is just one of the landshaping processes on Earth. Other mechanisms, 
such as mechanical erosion, would tend to increase the total denudation rate of a surface 
{Fleurant et al, 2008). The dissolved solids can take part in the crystallization of soluble 
deposits at the surface (e.g. calcretes or salt crusts in arid environments) or in the sub¬ 
surface (e.g. cements), when they saturate the liquids. They do not necessarily crystallize 
at the location where they dissolved {Ford and Williams, 2007). 

4 Denudation rates to determine ages of terrestrial karstic 
landscapes 

Using the denudation rates, theoretical timescales for the development of dissolution land- 
forms of a given depth can be inferred. We test hereafter this hypothesis on terrestrial 
instances by comparing ages determined from the denudation rates and ages determined 
from relative or absolute chronology. For carbonates, the partial pressure of CO 2 is 
computed based on the MOD16 Global Evapotranspiration Product of the University 
of Montana/ESRI Mapping Center. The notation GEyrs, MEyrs and kEyrs will refer to 
timescales expressed respectively in Giga Earth years (10® years). Million Earth years (10® 
years) and kilo Earth years (10^ years) in the following sections. 

4.1 Example under a hyper-arid climate: the caves of the Mount Sedom 
diapir (Israel) 

Mount Sedom is a 11 x 1.5 km salt diapir located in Israel, 250 meters above the level 
of the Dead sea, under a hyper-arid climate {Frumkin, 1994, 1996). It is composed of a 
Pliocene-Pleistocene (< 7 MEyrs) halite basis covered by a 5 - 50 m-thick layer composed 
of anhydrite and sandstones. Several sinkholes and caves exist in Mount Sedom {Frumkin, 
1994). 

In this part of Israel, the present-day hyper-arid climate leads to average rainfall rates 
of 50 mm/Eyr, of which 10 - 15 mm/Eyr constitutes the average effective rainfall rate 
(r). Gvishim and Mishquafaim Caves are respectively 60 m and 20 m-deep caves and are 
essentially carved into halite in the northern part of Mount Sedom. Radioisotopic age 
measurements in several caves of Mount Sedom give an Holocene age, about 8 kEyrs (and 
around 3.2 - 3.4 kEyrs from age measurements performed in Mishquafaim Cave) {Frumkin, 
1996). 

Assuming that underground dissolution is somewhat connected to the rainfall percolat¬ 
ing through anhydrite and reaching halite, the denudation rates computed with our model 
(Equation 1) over this region are DR ~ 1.69 — 2.53 mm/Eyr. The timescales required to 
form these caves would be between 23.7 and 35.5 kEyrs for Gvishim Cave, and between 
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7.9 and 11.8 kEyrs for Mishquafaim Cave, longer than those estimated from radioisotopic 
measurements. 

However, it should be noted that Mount Sedom experiences a rather deep dissolution, 
the load of dissolved solids in the underground liquid water flowing in its caves being up to 

3 times that of surface runoff {Frumkin, 1994). Therefore, considering that the denudation 
rate should be 3 times greater, we find timescales of development between 7.9 and 11.8 
kEyrs for Gvishim Cave and between 2.6 and 3.9 kEyrs for Mishquafaim Cave, in better 
agreement with age measurements. 

4.2 Example under a semi-arid climate: The Etosha super-pan (Namibia) 

The Etosha Pan, located in Namibia under a semi-arid climate, is a flat 120x60 km karsto- 
evaporitic depression that has already been suggested as a potential analogue for Titan’s 
lacustrine depressions {Bourgeois et ai, 2008; Cornet et ai, 2012). The depression is about 
15 to 20 m-deep at most, and has been carved into the carbonate layer essentially composed 
of calcretes and dolocretes lying ontop a middle Tertiary-Quaternary detritical sedimentary 
sequence that covers the Owambo basin (Buck and Rose, 1996; Hipondoka, 2005). This 
sedimentary sequence is believed to have accumulated under a semi-arid climate, relatively 
similar to the current one. The age of the calcrete layer is not known precisely, but its 
formation is believed to have started during the Miocene/Pliocene transition about 7 
MEyrs ago {Buck, 1997; Buck and Trippner, 1997), or even later during the late Pliocene, 

4 MEyrs ago [Miller et al, 2010). The Etosha Pan is believed to have developed at the 
expense of the calcrete layer since 2 MEyrs [Miller et al., 2010). 

In the Owambo sedimentary basin, AET ~ 350 mm/Eyr, yielding Pco 2 — 3 matm. 
Precipitation rates during the summer rainy season reach up to 500 mm/Eyr, which leads 
to r = 150 mm/Eyr at most. These conditions lead to DR = 6.5 /rm/Eyr in a substra¬ 
tum composed of calcite (calcrete) and DR = 5.6 /rm/Eyr in a substratum composed of 
dolomite (dolocrete). These denudation rates give an approximate age for the Etosha Pan 
of about 2.31 - 3.08 MEyrs in calcretes and 2.68 to 3.57 MEyrs in dolocretes. These ages 
are consistent with ages estimated from geological observations (< 4 MEyrs). 

4.3 Example under a temperate climate: the doline of Crveno Jezero 
(Croatia) 

Crveno Jezero is a 350 m-wide collapse doline located in the Dinaric Karst of Croatia, 
under a presently mediterranean temperate climate. This part of the Dinaric Karst is 
essentially composed of limestones and dolomites that have been deposited during the 
Mesozoic (Triassic to Cretaceous), when the area was covered by shallow marine carbonate 
shelves [Vlahovic et al., 2002; Mihevc et al., 2010). The cliffs forming the doline are up 
to 250 m in height. The bottom of the depression is covered by a lake, named Red Lake, 
about 250 m-deep. The total depth of the doline is therefore about 500 m [Garasic, 2001; 
Garasic, 2012). 

Although the formation process of the collapse dolines in the area is still not completely 
understood, especially in terms of the relative importance and timing of the collapse 
compared to the dissolution in the development of the structure, the uplift linked to the 
formation of the Alps (Eocene-Oligocene) led to the exposure of these carbonates and 
their subsequent karstification since at least Oligo-Miocene times (last 30 MEyrs) [Sket, 
2012). At some point, dissolution occurred to form an underground cave; then the load 
of the capping rocks exceeded their cohesion, leading to the brutal or progressive collapse 
of the cap in the empty space beneath. The timescale calculated for the doline formation 
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therefore constitutes a higher estimate. Age determinations in the northern Dinaric caves 
have been performed and show that karstification would be older than 4 to 5 MEyrs {Zupan 
Hajna, 2012). Crveno Jezero thus developed between 4 and 30 MEyrs ago. 

In this region of Croatia, AET ~ 500 mm/Eyr, which leads to Pco 2 — matm. 
Annual mean precipitation rates are about 1100 mm/Eyr (Mihevc et al, 2010), leading to 
T = 600 mm/Eyr. These parameters give DR = 31.8 /rm/Eyr in limestones dominated by 
calcite {DR = 27.5 /rm/Eyr in limestones dominated by dolomite). The time required to 
dissolve 500 m of calcite under these present conditions would be around 15.7 MEyrs (18.2 
MEyrs in dolomite). It is likely that cap rocks fell into the cave, leading to a younger age 
(1 m of calcite and dolomite dissolve in about 31 and 36 kEyrs respectively). However, 
the age determined by our method is still consistent with the geological records. 

4.4 Example under a tropical climate: the Xiaozhai tiankeng (China) 

The Xiaozhai tiankeng is a 600 m-wide collapse doline located in the Chongqing province of 
China, under a tropical climate. The depth of this doline is evaluated between 511 and 662 
m depending on the location (Xuewen and Weihai, 2006; Ford and Williams, 2007). It is 
organized into two major collapse structures, an upper 320 m-deep structure and a lower 
342 m-deep shaft. The tiankeng developed into Triassic limestones in a karst drainage 
basin. It results from various surface and subsurface processes such as dissolution, fluid 
flows and collapse. Tiankengs of China would have formed during the late Pleistocene 
(last 128,000 Eyrs) {Xuewen and Weihai, 2006) under an equatorial monsoonal climate 
similar to nowadays, emplaced 14 MEyrs ago {Wang and Li, 2009). 

In the area, annual precipitation rates are about 1500 mm/Eyr, for an annual evapo- 
transpiration rate of about 800 mm/Eyr, which leads to Pc 02 — 0.0123 atm and r = 700 
mm/Eyr. Under these conditions, DR = 49.6 ^um/Eyr (calcite) or DR = 42.6 /um/Eyr 
(dolomite) and the timescale to form the Xiaozhai tiankeng would be between 10.3 and 
13.3 MEyrs (calcite) and 12.0 and 15.5 MEyrs (dolomite), considerably longer timescales 
than those estimated from geological records. However, such a difference is to be expected, 
since the formation of gigantic tiankengs are subject to complex interconnected surface 
processes not taken into account in these simple calculations (collapse and underground 
circulation of water). This illustrates the limits of our method, generally over-estimating 
the ages of karstic features when they do not result from dissolution only. 

5 Denudation rates on Titan 

Equation 1 has to be adapted for Titan. Erom thermodynamics, we compute molar vol¬ 
umes (Um,i in m^/mol, see Appendix C) and mole fractions at saturation (Xj^satj see 
Appendix B) instead of mass densities and molalities. We considered the Ideal and the 
non-ideal Regular Solutions Theories (1ST and RST respectively) and ensured that our 
results are consistent with experimentally determined solubilities. The RST has been 
developed to determine the solubility of non-polar to slightly polar molecules, such as 
simple hydrocarbons or carbon dioxide, in non-polar solvents. It is less likely to be ap¬ 
propriate for polar molecules such as nitriles/tholins and water ice. We therefore always 
provide a comparison between the 1ST and the RST in order to assess the uncertainty 
of the calculation by using the RST for polar molecules (see also Appendix B). However, 
we consider the RST, which provides our lowest and more realistic solubilities (at least 
for hydrocarbons), as the most suited model for the study. We assume that dissolution 
is instantaneous, which is not unreasonable given the rapid saturation of solid hydrocar- 
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bons in liquid ethane inferred from recent dissolution experiments [Malaska and Hodyss, 
2014), compared to the long geological timescales considered in our work. For a binary 
solute-solvent system, Equation 1 becomes: 

X- + ■ 

j-^ T-t ^^ 2 ,sat m ,2 

= -T7- —f -5 

1 - Xi sat , 

m,solv 

where and are the molar volumes of the solid (crystallized) phase and of 

the liquid phase respectively. 

5.1 Comparison between Titan and the Earth over terrestrial timescales 

In order to compare the behavior of Titan’s solids in Titan’s liquids with that of min¬ 
erals in liquid water on Earth, we computed Titan’s and Earth’s denudation rates over 
Earth timescales by varying the precipitation rate between 0 and 2 m/Eyr. This range 
of precipitation rates encompasses the expected range on Titan (up to 1.2 - 1.3 m/Eyr 
during the rainy season {Schneider et al, 2012)). It also covers the precipitation range 
between arid and tropical climates on Earth. This comparison is illustrated in Eigures 
3 (for hydrocarbons and carbon dioxide ice) and 4 (for nitriles and water ice), for which 
we fixed Titan’s temperature to 91.5 K (the surface temperature during the rainy season 
according to Schneider et al. (2012)) and Earth’s temperature to 298.15 K (with a partial 
pressure of carbon dioxide equal to 0.33 matm). 

All organic compounds except CnHuN would behave like common mineral salts 
(halite, gypsum or anhydrite) according to the 1ST, which means that they would experi¬ 
ence dissolution rates in the range of hundred //m to a few cm over one Eyr. According to 
the RST and assuming that the liquid is methane, C 2 H 4 would behave like halite (cm-scale 
dissolution), C 4 H 10 and C 4 H 4 like gypsum (hundred fim to mm-scale dissolution), C 2 H 2 , 
C 3 H 4 , C 4 H 6 and CO 2 like carbonates (up to hundred /rm-scale dissolution). CgHg and 
nitriles would be less soluble in methane than calcite is in pure liquid water, with the least 
soluble nitriles being CnHuN, HCN and CH 3 CN. As expected, water ice is completely 
insoluble according to the RST, developed for non-polar molecules. 

The denudation rates of all organic compounds are higher in ethane and propane than 
in methane. Nitriles and CeHg, which are poorly soluble in methane, are quite soluble in 
ethane and propane (except CnHuN, HCN and CH 3 CN still under the calcite level) and 
reach denudation rates similar to those of carbonates or even gypsum. C 2 H 2 and C 3 H 4 
would behave like carbonates (in ethane) or gypsum (in propane). C 4 H 4 and C 4 H 6 behave 
like salts in those liquids and are even more soluble than gypsum. C 4 H 10 and C 2 H 4 are 
halite-like materials, extremely soluble in ethane and propane. Therefore, the likelihood of 
developing dissolution landforms in a hydrocarbon dominated substrate, by analogy with 
the Earth, is high. 

5.2 Present denudation rates on Titan over a Titan year 
5.2.1 Net precipitation rates on Titan 

Computing denudation rates over Titan timescales requires us to define the evolution of 
the net precipitation over one Titan year (1 Tyr ~ 29.5 Earth years). Titan’s climate is 
primarily defined by a rainy warm “summer” season and a cold dry “winter” season, both 
spanning about 10 Earth years. Southern summers are shorter and more intense than 
those in the north {Aharonson et al, 2009). Precipitation occurs as sporadic and intense 


[m/Eyr] 


(3) 
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rainstorms during summer, when cloud formation is observed {Roe et al, 2002; Schaller 
et ai, 2006, 2009; Rodriguez et al., 2009, 2011; Turtle et al., 2011a,b). 

A few Global Circulation Models (GCM) attempt to describe, at least qualitatively, 
the methane cycle on Titan (e.g. Rannou et al. (2006), Mitchell (2008), Tokano (2009) and 
Schneider et al. (2012)). Usually, net accumulation of rain is predicted at high latitudes > 
60° during summer, in agreement with the presence of lakes, whereas mid-to-low latitudes 
experience net evaporation, in agreement with the absence of lakes and the presence of 
deserts. Quantitatively, the model predictions are subject to debate since they depend 
on their physics (e.g. cloud microphysics, size of the methane reservoir, radiative transfer 
scheme). Still, they remain our best estimates about Titan’s climate. 

During the summer season, the model of Rannou et al. (2006) predicts net precipitation 
rates of methane lower than 1 cm/Eyr at 70° latitude and up to 1 m/Eyr poleward of 
70°, equivalent of a few to 2.7 mm/Earth day (Eday) respectively. The model of 
Mitchell (2008) predicts precipitation rates of about 2 mm/Earth day (Eday) in the polar 
regions and along an “Inter-Tropical Convergence Zone” (ITCZ), nearly moving “pole-to- 
pole”. The intermediate and moist models of Mitchell et al. (2009) predict precipitation 
roughly varying between 2 and 4 mm/Eday at the poles and along the ITCZ. These rates 
are consistent with those estimated in Mitchell et al. (2011) in order to reproduce the 
tropical storms seen in 2010 {Turtle et al, 2011a) and with those estimated by the model 
of Schneider et al. (2012). The model of Tokano (2009) also predicts similar precipitation 
rates (800 to 1600 kg/m^ in half a Titan year, equivalent to precipitations between 3 and 
6 mm/Eday). 

Here, we use the methane net precipitation rates extracted from the GCM of Schneider 
et al. (2012) to compute the present-day denudation rate on Titan. Figure 5 represents 
the mean net precipitation rates at various polar latitudes. High latitudes > 80° would 
be quite humid (r = 7 — 8 m/Tyr). Lower latitudes would be less humid (r = 3 — 3.6 
m/Tyr at 70°, decreasing to r = 0.4 — 1.6 m/Tyr at 60°). Southern low latitudes would 
be much drier than northern latitudes over a Titan year as a result of sparser but more 
intense rainstorms. 

5.2.2 Case of pure compounds 

Figure 6 illustrates the denudation rates of a surface composed of pure organic compounds 
exposed to methane rains at several southern and northern polar latitudes according to 
the 1ST and RST hypotheses. Ethane and propane are not shown since the model of 
Schneider et al. (2012) only considers methane, but the behavior of Titan’s solids in these 
liquids is already discussed in Section 5.1. 

Over one Titan year, the denudation rates are the highest at high latitudes and the 
lowest at low southern latitudes. According to the 1ST, all compounds would experience 
dissolution on the order of a few mm to a few meters at almost all latitudes (salt-like 
material), except CnHuN, the denudation rate of which would be a few hundred nm. 
The dissolution of C 2 H 8 N 2 , CO 2 , HCN and CeHg at 60°S however would be on the order 
of several hundred /xm (carbonate-like material). 

According to the RST, C 2 H 4 and C 4 H 10 denudation rates are between a few mm 
up to a few meters per Titan year (salt-like materials). The dissolution of C 4 H 4 , C 4 H 6 , 
C 2 H 2 , CO 2 , C 3 H 4 is /xm to hundred /xm-scale over a Titan year (carbonate-like materials). 
The denudation rates of nitriles and CeHg are between 10“'^ — 10“^® mm/Tyr (for HCN, 
CH 3 CN, CiiHiiN) and a few ^um/Tyr (carbonate to siliceous-like materials). 

At high latitudes > 70°, we do not see much differences between the northern and the 
southern denudation rates, as expected from the similarities in precipitation rate between 
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the two poles shown in Figure 5. At low southern latitudes, net precipitation rates are too 
low over one Titan year to allow a rapid and significant dissolution. Interestingly, Ontario 
Lacus and Sikun Labyrinthus, two landforms compared with terrestrial karsto-evaporitic 
and karstic landforms (Malaska et al, 2010; Cornet et al, 2012), are observed at latitudes 
greater than 70°S, and no other well-developed dissolution-related landforms are seen at 
lower southern latitudes. 

Therefore, if Titan’s surface is composed of pure hydrocarbons, dissolution processes 
are likely to occur but the formation of a karstic-like landscapes would be roughly 30 times 
slower on Titan than on Earth due to Titan’s seasonality in precipitation. Of course, this 
latter consideration depends on the actual composition of the surface, which is unlikely to 
be pure, and of the accuracy of the climate model used. 

5.2.3 Case of a mixed surface layer 

We now assume the presence of a surface layer, the composition of which is proportional to 
the accumulation rates at the surface (hi) of solids coming from the atmosphere, calculated 
in the same way as Malaska et al. (2011) did: 

hi = Pi [m/Eyr] (4) 

where pi is the production rate of molecules (in molecules/m^/Eyr) listed in Table 1, 
is the molar volume of the solid (or subcooled liquid if the former is not known, in 
m^/mol) and Na is the Avogadro number (~ 6.022 x 10^^ mol“^). The composition of the 
mixed organic layer is then determined as percentages (/*) of each organic compound in 
the layer (so that fi = hi/ Y/i hi). This method allows us to consider the volume occupied 
by each molecule, which is important especially for tholins because these contribute up to 
20 % on average to the total production rates of molecules, but they build up to ~ 50 % 
of the total thickness of the surface deposits due to their higher molar volumes compared 
to those of simple hydrocarbons or nitriles. 

We consider 3 cases for the surface composition: without tholins and with C 2 H 8 N 2 or 
CiiHiiN as tholins. Tholins production rate is that mentioned in Cabane et al. (1992). 
Methane precipitation rates are those of Schneider et al. (2012). The denudation rates for 
these mixed organic layers (DRmix) are computed in a linear mixing model scheme where: 

D Rmix — E f^DR,. [m/Eyr] (5) 

i 

Figure 7 gives the repartition of denudation as a function of latitude and photochemical 
models at the end of a Titan year. According to the 1ST, the denudation rate of all mixed 
layers would be on the order of a few cm to a few meters over a Titan year (salt-like layers). 
According to the RST, the organic layer originating from the Lavvas et al. (2008b) model 
would be the most soluble (dissolution rates of a few cm to a few dm per Titan year, 
salt-like layer). The organic layers originating from the other models would be more 
carbonate-like layers over a Titan year (dissolution rates of a few tens of pm to a few mm 
per Titan year), whether tholins are included or not. The lowest solubility of all mixed 
layers is reached using a Krasnopolsky (2009)-type composition. Over a Titan year, the 
likelihood of developing dissolution-related landforms is therefore non-negligible, even if 
the surface is not composed of pure soluble simple solids. These mixed organic layers 
would behave like carbonate or salty terrestrial layers over a Titan year. 
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6 Discussion: How old are Titan’s karstic landscapes? 


Despite the strong assumptions of the method described in Section 4 to infer timescales 
of formation of terrestrial karstic landforms (we consider only chemical erosion at equi¬ 
librium without significant climate changes over the past few MEyrs), the resulting ages 
are consistent with ages determined by relative or absolute chronology or constitute upper 
limits. Therefore, the determination of the age of the lacustrine depressions on Titan 
probably result in maximum timescales of development. Titan’s climate is believed to 
have remained quite stable over the recent past, with a small periodic insolation variation 
of ±2 W/m^ at the poles during the last MEyr, for the current low insolation at the North 
pole [Aharonson et al, 2009). This probably brings some stability to the calculations. 

By applying our simple model, we compute the timescales needed to form a 100 m-deep 
depression by dissolution of a superhcial mixed organic layer under the current climate 
conditions evaluated by Schneider et al. (2012). These are shown in Table 3 and Eigure 
8 . We compute our formation timescales using both the 1ST and the RST since the 
solubility of polar molecules is not well constrained by the RST and could get closer to 
that computed using the 1ST. However, we consider timescales evaluated using the RST as 
our references since they present the most conservative values for the age of the lacustrine 
depressions. 

Independent of the thermodynamic theory considered. Titan’s lacustrine depressions 
would be young. Among all the compositions tested, the time needed to carve a 100 
m-deep depression by dissolution under current climate conditions at latitudes poleward 
of 70° would be between a few kEyrs (1ST) and 56 MEyrs (RST). At 60°N, a 100 m- 
deep depression would be created in 7.7 kEyrs (1ST) to 104.6 MEyrs (RST) while the 
same depression would be created in 27.6 kEyrs (1ST) to 375.1 MEyrs (RST) at 60°S. 
This strong difference between the two hemispheres could explain why Titan’s south polar 
regions are deprived of well developed lacustrine depressions compared to the North. 

It should be noted that the hypothesized timescale difference between the northern 
and the southern low latitudes results from an extrapolation of Titan’s current climate to 
the past. Aharonson et al. (2009) showed that Croll-Milankovich-like cycle with periods 
of 45 (and 270) kEyrs could exist on Titan, resulting in a N-S reversal in insolation 
and likely subsequent climate conditions. Over geological timescales, the N-S differences 
in denudation rates and timescales estimated at these latitudes could be smoothed by 
these cycles. However, as noted earlier, we made the same extrapolation for the Earth, 
whose climate dramatically changed over time, without obtaining unreasonable formation 
timescales. 

In any case, all these timescales are consistent with the youth of Titan’s surface as 
determined from: 

1. crater counting (0.3 — 1.2 GEyrs, Neish and Lorenz (2012)), 

2. dune sediment inventory (50 — 730 MEyrs, Sotin et al. (2012) and Rodriguez et al. 
(2014)), 

3. the flattening of the poles due to the substitution of methane by ethane in clathrates 
(500 MEyrs if restricted to the poles or 0.3 — 1.7 GEyrs if not, Choukroun and Sotin 
( 2012 )) 

4. the possible methane outgassing event (1.7 — 2.7 GEyrs, Tobie et al. (2006)). 

In summary, the morphology of Titan’s lacustrine depressions suggests that dissolution 
occurs on Titan. The denudation rates of pure organic compounds and a mixed organic 
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layer as compared to those of soluble minerals on Earth also supports this hypothesis. The 
timescales needed to dissolve various amounts of material as compared to the timescales 
of development of karstic landforms on Earth are also quite consistent in the sense that 
karstic landscapes are usually relatively young landscapes. Finally, the latitudinal repar¬ 
tition of denudation rates and timescales of dissolution is consistent with the latitudinal 
repartition of the possible dissolution-related landforms at the surface of Titan. The sur¬ 
face dissolution scenario for the origin of Titan’s lakes appears very likely and Titan’s 
lakes could be among the youngest features of the moon. 

7 Conclusion 

Titan’s lakes result from the filling of topographic depressions by surface or subsurface 
liquids. Their morphology led to analogies with terrestrial landforms of various origins 
(volcanic, thermokarstic, karstic, evaporitic or karsto-evaporitic). The karstic/karsto- 
evaporitic dissolution scenario seems to be the most relevant, given the nature of surface 
materials on Titan and its climate. We constrained the timescales needed for the forma¬ 
tion of Titan’s depressions by dissolution, on the basis of the current knowledge on the 
development of terrestrial karsts. 

We computed solutional denudation rates from the theory developed by White (1984). 
This simple theory needs three parameters: the solubilities and the densities of solids and 
liquids at a given temperature and a climatic parameter linked to the net precipitation rates 
onto the surface. We computed the solubilities of terrestrial minerals in liquid water at 
25° C and tested the model by computing the denudation rates and timescales of formation 
of several terrestrial examples of karstic landforms. 

We then applied the same model to Titan. We computed the solubilities of Titan’s 
surface organic compounds in pure liquid methane, ethane and propane at 91.5 K using 
different thermodynamic theories. We evaluated the molar volumes of liquid and solid Ti¬ 
tan’s surface compounds at 91.5 K and we used the results of the recent GCM of Schneider 
et al. (2012) as input for the precipitation rates of methane on Titan, which allowed us 
to compute denudation rates at several latitudes. Denudation rates have then been com¬ 
puted for pure organic compounds at Earth and Titan timescales and have been compared 
to those determined for soluble minerals on Earth. We also computed denudation rates 
for three different compositions of the surface organic layer. Over one Titan year, these 
mixed layers of organic compounds behave like terrestrial salts or carbonates, which in¬ 
dicates their high susceptibility to dissolution, though these processes would be 30 times 
slower on Titan than on Earth due to the seasonality (rainfall occurs only during Titan’s 
summer). 

We computed theoretical timescales for the formation of 100 m-deep depressions in 
mixed organic layers under present climatic conditions. As with dissolution landforms 
on Earth, Titan’s depressions would be young. At high polar latitudes, we found that 
the timescales of development for depressions are relatively short (on the order of 50 
MEyrs at maximum to carve 100 m) and consistent with the young age of Titan’s surface. 
These timescales are consistent with the existence of numerous lacustrine depressions and 
dissected landscapes at these latitudes. At southern low latitudes, the computed timescales 
are as long as 375 MEyrs due to the low precipitation rates. This low propensity to develop 
depressions by dissolution is consistent with their relative absence/bare formation at low 
latitudes. Over geological timescales greater than those of Titan’s Croll-Milankovitch 
cycles (45 and 270 kEyrs), this difference would probably be strongly attenuated. However, 
climate model predictions are not presently available over geologic timescales, and the 
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present-day seasonal climate variations are the best that can be currently constrained. 

The results of these simple calculations are consistent with the hypothesis that Titan’s 
depressions most likely originate from surface dissolution. Theoretical timescales for the 
formation of these landforms are consistent with the other age estimates of Titan’s surface. 
Future works could include the effects of rain in equilibrium with the nitrogen, ethane and 
propane atmospheric gases in the raindrop composition (e.g. Graves et al. (2008) or Glein 
and Shock (2013)). Experimental constraints on the solubility of gases and solids in liquids 
thanks to recent technical developments for Titan experiments {Luspay-Kuti et al., 2012, 
2014; Malaska and Hodyss, 2014; Chevrier et al, 2014; Leitner et al, 2014; Singh et al., 
2014) would also be of extreme importance for such work. Finally, the influence of other 
landshaping mechanisms such as collapse or subsurface fluid flows, which play a significant 
role in the development of some karstic landforms on Earth, could also be implemented. 
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A List of compound names for Titan 

CH 4 : methane (liquid) 

C 2 H 6 : ethane (liquid) 

CsHg: propane (liquid) 

C 4 H 10 : n-butane 
C 2 H 2 : acetylene 
C 2 H 4 : ethylene 
C 3 H 4 : methyl-acetylene 
C 4 H 4 : vinyl-acetylene 
C 4 H 6 : 1-3 butadiene 
CgHg: benzene 
HCN: hydrogen cyanide 
C 2 N 2 : cyanogen 
CH 3 CN: acetonitrile 
C 2 H 3 CN; acrylonitrile 
C 2 H 5 CN; propionitrile 

C 2 H 8 N 2 : 1-1 dimethyl-hydrazine (tholins-like) 

CiiHiiN: quinoline (tholins-like) 

H 2 O: water ice Ih 
CO 2 : carbon dioxide ice 


14 


B Solubility of solids on Earth and Titan 

B.l Dissolution processes on Earth 

Dissolution occurs in karstic to evaporitic areas on Earth, where the dominant minerals 
are calcite, dolomite, halite, anhydrite or gypsum. Basic aqueous reactions of congruent 
dissolution for these common soluble minerals (i.e., all components of the minerals dissolve 
entirely) are summarized below: 

CaCOs + H 2 O+ CO 2 ^ Ca2+ + 2 HCO^ 

CaMg(C03)2 + 2 H 2 O + 2 CO 2 ^ Cd?+ + Mg2+ + 4 HCO 3 
CaS 04 . 2 H 20 ^ Ca2+ + SO^” + 2 H 2 O 
CaS 04 + H 2 O ^ Ca2+ + SO^” + H 2 O 
NaCl+ H 2 O ^ Na+ + Cl" + H 2 O 

Halite (NaCl), gypsum (CaS 04 . 2 H 20 ) and anhydrite (CaS 04 ) dissolve by pure dis¬ 
sociation, and calcite (CaCOs) and dolomite (CaMg(C 03 ) 2 ) dissolve by dissociation and 
acid dissolution, i.e. with the help of the CO 2 gas dissolved in water {Langmuir, 1997; 
Ford and Williams, 2007; Brezonik and Arnold, 2011). 

The physical and thermodynamic parameters needed to compute solubilities are given 
in Table 4. Ideal and non-ideal electrolyte theories are considered. We also consider the 
acid dissolution of carbonate minerals. 

B.1.1 Pure dissociation of minerals: halite, anhydrite and gypsum 

For a given chemical reaction implying A and B and producing C and D, with their 
respective stoichiometric numbers a, h, c and d, the law of mass action gives: 

aA + bB <;=> cC + dD, (6) 

where one can define a thermodynamical equilibrium constant (or thermodynamic 
solubility product) for the reaction K^q (often written pKeq, with pKeq = —^ogKgq), so 
that: 

^ {CriD)^^ ^ {mcTcrimDrD)<^ 

(H)“(H)^ {mAT a)'" (mBT 

where (i) denotes the activity of the ith compound, being the product of an activity 
coefficient Tj and a molality ruj (in mol/kg of solvent). For liquids and solids reacting 
together (or being produced) during a congruent dissolution reaction, the activity is set to 
unity. The equilibrium constant can be derived from thermodynamics at standard state 
(T = 25°C and P = 1 atm) using the Gibbs-Helmholtz equation: 

ArG° =-RT InK^q =-2.303 RT logKeq, [J/mol] ( 8 ) 

where Aj-G° is the standard Gibbs free energy of reaction (in J/mol), calculated using 
the Hess’ law: 


A,G°= niAfG°- [J/mol] (9) 

products reactants 

with AfG° being the standard Gibbs free energy of formation of each species (in J/mol) 
and n their stoichiometric numbers. Following the same principle, we compute the stan¬ 
dard enthalpy of reaction A^H° (J/mol) from the individual enthalpies of formation Af77° 
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(in J/mol). The standard enthalpy of reaction is then used to compute the equilibrium 
constants at temperatures that differ from standard conditions using the Van’t Hoff equa¬ 
tion {Langmuir, 1997) as follows: 


In 


Kt2 

Kti 


R T 2 J ’ 


( 10 ) 


where Kti and Kt 2 are the equilibrium constant of a given chemical reaction at two 
different temperatures Ti and T 2 (in K) (the subscript 1 is used to indicate the reference 
state, 25°C in this case). Table 5 gathers the thermodynamic parameters for all the 
considered dissolution reactions. It also displays each equilibrium constant equation for 
common terrestrial minerals. 


B.1.2 Acid dissolution of carbonates: calcite and dolomite 

For calcium and magnesium carbonates, a set of reactions happens, not only involving 
the mineral dissociation itself {Kd or Km), but also including the carbon dioxide gas 
dissolution in water {Kco 2 )j which produces carbonic acid that acidifies water. Then, 
carbonic acid rapidly dissociates into bicarbonate ions (Ki), which are also created by 
the association of protons H'*' and carbonate ions ( 7 ^ 2 ). The thermodynamic properties 
of these intermediate reactions to the dissolution of calcite and dolomite are summarized 
in Table 5. The activity of carbon dioxide is approximately equal to its partial pressure 
{Pc 02 , see Langmuir (1997) and Ford and Williams (2007)) so that, for calcite acid 
dissolution, one can write: 


_ (Ca^+) {HCO^)^ _ Ki Kco2 K^ 
PCO2 K2 

and for dolomite acid dissolution: 

_ (Ca2+) {Mg^+) (HCO^)^ _ Kf K^^, Km 

PU 


( 11 ) 


( 12 ) 


Assuming a given partial pressure of carbon dioxide, one can compute the activity, 
molality and mole fraction at saturation of calcite and dolomite in water. The molality of 
CO 2 in the system is calculated thanks to the Henry law {mco 2 — kniP) PCO 2 , Kh{T) 
being Henry’s law constant, varying with temperature). We performed the calculations 
for 3 values of Pc 02 25°C: Pc 02 — 0-33 matm, which represents the normal dry air 
{Langmuir, 1997) that could be quite analogous to atmospheric conditions in arid/semi- 
arid areas, and Pco 2 — 0.01—0.11 atm, which represents values encountered in more humid 
areas such as under equatorial/tropical climates {Ford and Williams, 2007; Fleurant et ai, 
2008). 


B.1.3 Activity coefficients for electrolyte solutions 

We infer the molalities m* and activities (i) of ions in solution from Keq at different 
temperatures under an Ideal (1ST, Tj = 1) and a non-ideal Electrolyte (EST, T* / 1) 
Solution Theory. Eor the EST, the activity coefficients Tj are calculated by iteration using 
the extended Debye-Hiickel equation as modified by Truesdell and Jones (1974): 


log Ti 


-A; 


VI 


1 -j- p ^/l 


+ bid, 


(13) 
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where I is the molal ionic strength of the solution, defined as follows: 

1 ” 

^ = 2 YI [mol/kg] (14) 

i=l 

Zi is the charge of i, ai and bi are the 2 parameters of the Truesdell and Jones (1974) 
equation (a modihed hydrated radius, in A, and a purely empirical parameter respectively). 
A and B are two variables depending on the temperature T (in °C), defined as follows 
{Ford and Williams, 2007): 


A = 0.4883 + 8.074 X 10"^ T, (15) 

B = 0.3241 + 1.6 X lO-'^T. (16) 

Results from these calculations are given in Table 6. The calculation of activity co¬ 
efficients has a relatively minor impact in the change of molality for carbonates, while a 
clear difference can be seen for salts (anhydrite and gypsum). The presence of CO 2 in the 
carbonate-water system greatly increases the amount of dissolved carbonates, until reach¬ 
ing solubilities roughly similar to those of gypsum and anhydrite. The EST molalities has 
been used for all minerals except for halite, which possesses a ionic strength too strong to 
be computed using our Debye-Hiickel equation. 


B.2 Solubility of Titan’s solids 

The solubility of Titan’s solids in pure liquid hydrocarbons is calculated according to the 
Van’t Hoff equation: 

InTjXj^sat — It ■ T/ ’ 

where T, is the activity coefficient of the solute, its mole fraction at saturation, 

AinH° its enthalpy of melting (in J/mol) evaluated at Tm,i, its melting temperature (in 
K), which differs from the temperature of the solution T (K). R is the ideal gas constant. 
Melting properties of Titan compounds are given in Table 7. 

We consider the case of the Ideal Solutions Theory (1ST), for which Xi^sat is calculated 
by assuming Tj = 1 (all the molecules share the same affinity with each other), and the 
case of the Regular Solutions Theory (RST) of Preston and Prausnitz (1970) for which 
Tj A 1 (molecules have preferential affinity with one or the other). The latter approach 
has already been used in several publications dealing with the lakes and seas composition 
{Raulin, 1987; Dubouloz et ai, 1989; Cordier et al, 2009, 2013b). 

The RST requires to compute Tj’s. Assuming the subscripts 1 for the solvent and 2 
for the solute, the RST gives: 


InTa 


(dl — 62 )'^ + 2 

RT ”’2' 




(18) 


where 14^2 is the subcooled liquid molar volume of the solute (in m^/mol), I 12 is the 
empirical interaction parameter between the solvent and the solute taken from various 
sources (e.g. Preston and Prausnitz (1970) and Szczepaniec-Cieciak et al. (1978)), (^1 is 
the volume fraction of the liquid, computed as follows: 


4>i 


( 1 -^ 2 ) R ^,1 

(1-^2)R^,1+W2R^,2’ 


(19) 
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and 8 i and 82 are the Hildebrand solubility parameters of the solvent and the solute 
respectively (in computed from their enthalpy of vaporization A^Hi (J/mol), 

as follows: 

[(J/m=)‘/^l (20) 

^m,i 

where A^Hi is taken in the literature at the boiling point temperature Tf, (Table 7) for 
each compound and extrapolated down to lower temperatures using the Watson equation 
{Poling et ai, 2007): 

= A^H.^Tb (r^^) °( 21 ) 

The solubilities are given for a wide range of temperatures in Figures 9 (non-polar 
molecules) and 10 (polar molecules), ontop of which we also reported experimental data 
points gathered from the literature, htted by an empirical power law {R^ > 0.98) in order 
to extrapolate the experimental values at low temperatures for comparison. 

The 1ST hypothesis provides a hrst estimate of the solubility, as seen in Figures 9 and 
10. It does not consider the different behavior of the solutes in the different solvents, 
but presents the advantage of not depending on approximate calculations of the activity 
coefficients of each species. According to this theory, all simple hydrocarbons would thus 
be rather solubles at low temperatures, the least soluble being nitriles and CeHg. 

The RST hypothesis provides a second estimate and shows that solids are less soluble 
in methane than in ethane and propane. Whatever the liquid, simple nitriles generally 
show an increase in solubility with a higher number of C atoms in their carbon chain. 
C 2 H 4 , C 4 H 10 and C 2 H 2 are the most soluble compounds again, joined by C 4 H 4 , C 4 H 6 , 
CO 2 , C 3 H 4 and finally CgHg. Nitriles are the least soluble compounds. At cryogenic 
temperatures (< 100 K), the simulated solubilities are in good agreement with the power 
fits of experimental data (with the exception of n-butane, which shows a higher solubility 
in experiments). It is also worth noting that recent solubility experiments of benzene in 
liquid ethane have been performed at 94 K by Malaska and Hodyss (2014). The 
for benzene in ethane they determined experimentally (1.11 ± 0.09 x 10“®) is also quite 
close to our theoretical value computed using the RST at 94 K (4.42 x 10“®), though 
slightly lower. This difference is probably due to the dissolution of nitrogen into liquid 
ethane during the experiment, which tends to lower the solubility of benzene. Benzene 
also reached saturation so quickly in liquid ethane (in less than 2h) {Malaska and Hodyss, 
2014), compared to the timescales considered in the present work, that dissolution can be 
assumed instantaneous. Finally, as expected, the solubility of water ice is unconstrained 
given the considerable range between 1ST and RST estimates (about a factor 10^^), but 
is probably low. 

Solubilities calculated using the 1ST and RST are also given at 91.5 K in Table 8 along 
with those gathered from our empirical fits of experimental data and the literature. The fits 
of experimental data are quite consistent with our RST values. Previous studies document 
the solubility of pure solids in liqnid mixtures in equilibrinm with the atmosphere (thus 
composed of methane, ethane, propane and/or nitrogen) at various temperatures. Direct 
comparisons are therefore tentative, since the liquid composition, temperatures and the 
thermodynamic parameters sources are heterogeneous. Nitrogen tends to decrease the 
solnbility of solids in liqnids while ethane and propane tend to increase it. It should be 
noted that the solubility of hydrocarbons is quite in agreement between our study and 
others {Raulin, 1987; Dubouloz et ai, 1989; Glein and Shock, 2013) but the solubility of 
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nitriles appears lower in our simulations than in previous ones {Raulin, 1987; Dubouloz 
et ai, 1989; Cordier et ai, 2013a), for which they lie between our 1ST and RST estimates. 


C Molar volumes of solvents and solutes 

If the density of a compound at a given temperature (T) is known from experimental data, 
it is straightforward to derive its molar volume l^,i(T), as follows; 

Vra,i{T) = Mi/pi{T) [mVmol] (22) 

This is the case for water ice (Loerting et al, 2011) and minerals {Lide, 2010). However, 
experimentally determined densities for solids and liquids at the very low temperatures 
relevant to Titan are rather rare. We describe two complimentary techniques to evaluate 
the molar volumes of solids and liquids. 


C.l Subcooled liquid molar volumes 


We first evaluate molar volumes at a given temperature T using the Rackett equation 
{Spencer and Danner, 1972; Poling et al, 2007). This method is designed to estimate 
the saturated subcooled liquid molar volumes of pure hydrocarbons and organic solvents 



[mVmol] 


(23) 


R is the ideal gas constant (= 8.3144621 J/mol/K), Tc and Pc are the critical temper¬ 
ature (in K) and pressure (in Pa or J/m^) respectively, Zc is the critical compressibility 
factor, similar to the Rackett parameter (Zra = 0.29056 — 0.08775t<;, iv being the ac- 
centric factor, see Poling et al. (2007) for more details) and the reduced temperature 
(Tf^i = T/Tc^i). All these parameters are given in Table 9. An alternative form of this 
equation is given by: 

= Vc,i Zc,i [mVmol] (24) 


We investigated uncertainties in the molar volumes computed using the Rackett equa¬ 
tions (23 and 24), with the two different parameters (.^rA) taken from database or com¬ 
puted using the accentric factor, and Z^ computed from the ideal gas equation at the 
critical point). The results are reported in Table 10 and Figures 11 (liquids), 12 (non¬ 
polar solids) and 13 (polar solids). Differences in molar volumes are only significant for 
nitriles, going up to 2cj < 24 cm^/mol for propionitrile. To reduce possible over or under¬ 
estimates for nitriles molar volumes and solubilities, we use Equation 24 with Zra, which 
is our “intermediate” case and also the most accurate estimates of 14^^ {Poling et ai, 
2007). 


C.2 Titan’s solids molar volumes 

The solid molar volumes were determined from crystal structures available in the literature 
{Dulmage and Lipscomb, 1951; Dietrich et ai, 1975; Simon and Peters, 1980; Refson 
and Pawley, 1986; Antson et al., 1987; Etters and Kuchta, 1989; McMullan et al, 1992; 
Craven et ai, 1993). The CH 3 CN, HCN, C 4 H 10 and C 2 H 2 have a temperature dependent 
polymorphism. For each of these compounds, the selected crystal structure corresponds to 
the one observed at the temperature of Titan’s surface (i.e. 90 - 95 K). Except for CH 3 CN, 
the cell volume was experimentally determined for different temperatures. A second order 
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polynomial was used to fit the cell volume as a function of the temperature to extrapolate 
them at Titan’s temperatures (90 - 95 K). 

Since the pressure on Titan (around 1.5 bar) is not the same as on Earth (1 bar) where 
the crystal cell volume were determined, we checked the influence pressure on molar vol¬ 
umes by calculations based on the Density Functional Theory (DFT) using the VAST 5.3 
package. Cell volumes of each compound were optimized using vdW-DF2 {Lee et al, 2010; 
Klimes et al, 2011) functional involving a non-local kernel for the electronic correlation 
energy calculation. This functional allows to reproduce dispersion interactions (such as 
Van der Waals interactions) which are major inter-molecular interaction encountered in 
molecular crystals. The Fncut value, defining the basis set size, was fixed to 800 eV in 
order to suppress Pulay’s stress. The k-points sampling was done with a 12 x 12 x 12 
Monkhorst-Pack grid. The core electrons were described with the projector-augmented 
plane wave (PAW) approach. Calculations were performed with and without static pres¬ 
sure (set to 100 bar) to compare cell volumes. The larger variation was observed for C 4 H 10 , 
with a volume decrease of 0.44 % upon pressure. This is reasonably low to conclude that 
cell volumes measured on Earth (P=l bar) should be very similar to the ones on Titan 
(P=1.5 bar). 

We reported on Figures 12 and 13 the predicted crystalline molar volumes ontop of 
the experimental and Rackett values (individual values also given in Table 11). Overall, 
we get a quite good agreement between all these estimates. We estimate the solid molar 
volumes at 91.5 K using a simple empirical linear fit between experimental and predicted 
points. The Rackett method is used to give an approximate molar volume for the solids 
for which we do not have experimental data. 
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Figure 1: RADAR SAR mosaic of the TA (2004) to T95 (2013) flybys in north and south 
polar stereographic projections (down to 55° latitudes). All current lakes observed to date 
are located above 70° in latitude (blue circles). 
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Figure 2: RADAR SAR views of Titan’s second largest sea Ligeia Mare (a), empty 
and liquid-covered lacustrine depressions (b), and the heavily dissected terrains of Sikun 
Labyrinthus (c). Ligeia Mare clearly differs in size and shape from the smaller lacustrine 
depressions, which suggests a different formation mechanism. Lacustrine depressions and 
labyrinthic terrains pertain to the class of features potentially related to the dissolution 
of the surface. 
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Figure 3: Denudation rates (in mm/Eyr) of non to slightly polar molecules (hydrocarbons 
and carbon dioxide ice) in Titan’s liquids at 91.5 K (light curves: 1ST, bold curves: RST) 
compared to those of Earth’s minerals in liquid water at 298.15 K (EST). The partial 
pressure of carbon dioxide for acid dissolution is set to 0.33 matm. 
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Figure 4: Same as Figure 3 but for polar molecules (nitriles and water ice). 
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Figure 5: Net precipitation rates in mm/Earth day (Eday) at high southern and northern 
latitudes from Schneider et al. (2012) over a Titan year (29.5 Eyrs). SSS and NSS: 
Southern and Northern Summer Solstices. 


Ideal Solutions Theory Regular Solutions Theory 




Figure 6: Denudation rates (in mm/Tyr) of Titan’s surface solids in liquid methane accord¬ 
ing to the methane net precipitation rates of Sehneider et al. (2012) over an entire Titan 
year for the 1ST and RST approaches. In brown are represented simple hydrocarbons, in 
blue the nitriles and in green carbon dioxide. 
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Figure 7: Denudation rates (in mm/Tyr) of Titan’s surface solids in liquid methane ac¬ 
cording to the Schneider et al. (2012) model net precipitation rates over an entire Titan 
year at various latitudes and for different surface layer compositions established from pho¬ 
tochemical models, including (dashed lines) or not (solid lines) the influence of tholins 
(dashed line: with C 2 N 8 N 2 ; dashed-dot line: with CnHuN). 
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Figure 8: Same as Figure 7 but for the timescales of formation (in Eyr) of 100 m-deep 
depressions on Titan. 
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Figure 9: Saturation mole fraction of Titan’s “non-polar solids” (hydrocarbons and carbon 
dioxide ice) in Titan’s liquids computed using the RST and 1ST, and different Rackett 
equations (see Section C.l). The darkest curves represent Xi^sat evaluated by Equation 24 
with the Zra parameter. Diamonds are experimental values collected for C 4 H 10 {Brew, 
1977; Kuebler and McKinley, 1975), C 2 H 2 {Neumann and Mann, 1969), C 4 H 6 {Preston 
et ai, 1971), CgHg {von Szalghary, 1972; Luks et at, 1981; Kuebler and McKinley, 1995; 
Malaska and Hodyss, 2014) and CO 2 {Clark and Din, 1953; Davis et ai, 1962; Cheung 
and Zander, 1968; Preston et al, 1971) in the corresponding liquids. Dashed lines are 
empirical power fits of the experimental points (with > 0.98). 
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Figure 10: Same as Figure 9 but for Titan’s 


polar solids (nitriles and water ice). 
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Figure 11: Differences in the calculated molar volumes of Titan’s liquids due to the use of 
different Rackett equations and compressibility factors. 


Table 1: Production rates (in molecules/cm^/s) of simple solid organics in Titan’s atmo¬ 
sphere derived from several photochemical models. Not reported here is the production 
rate of tholins, estimated by Cabane et al. (1992), Rannou et al. (2003) and Krasnopolsky 
(2009) to be between 0.2 and 7 kg/cm^/GEyr (Giga Earth year). The list of compound 
names is available in Appendix A. Y84: Yung et al. (1984), T95: Toublanc et al. (1995), 
L96: Lara et al. (1996), W04: Wilson and Atreya (2004), LOS: Lavvas et al. (2008a,b), 
K09: Krasnopolsky (2009). 
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Figure 12: Same as Figure 11 but for “non-polar” solids (hydrocarbons and carbon dioxide 
ice). Red and blue diamonds are experimental and predicted values respectively. Sources: 
C 4 H 10 {Refson and Pawley, 1986), C 2 H 2 {McMullan et al, 1992), CeHg {Craven et ai, 
1993) and CO 2 {Simon and Peters, 1980). 
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Figure 13: Same as Figure 12 but for polar solids (nitriles and water ice). Sources: HCN 
[Dulmage and Lipscomb, 1951), CH 3 CN (Barrow, 1981; Antson et al, 1987) and H 2 O 
(Loerting et al, 2011). 
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Table 2: Solutional denudation rates in mm/Eyr for terrestrial minerals at 298.15 K 
and under various climatic conditions (CO 2 content and precipitation rates). All values 
correspond to the EST hypothesis except those for halite (1ST). 


Name 

T = 30 cm/Eyr 

r = 1 m/Eyr 

Dissolution by dissociation 

Halite 

50.643 

167.122 

Gyspum 

0.504 

1.662 

Anhydrite 

0.224 

0.740 

Calcite 

8.55x10""^ 

2.82x10"® 

Dolomite 

7.84x10""^ 

2.59x10"® 

Acid dissolution 

Calcite, PCO 2 = 0.33 matm 

0.006 

0.020 

Dolomite, PCO 2 = 0.33 matm 

0.005 

0.017 

Calcite, PCO 2 = 0.01 atm 

0.020 

0.066 

Dolomite, PCO 2 = 0.01 atm 

0.017 

0.057 

Calcite, PCO 2 = 0.11 atm 

0.047 

0.156 

Dolomite, PCO 2 = 0.11 atm 

0.040 

0.132 


Table 3: Timescales required to dissolve 100 meters of a surface organic layer (compo¬ 
sition given by photochemistry in the atmosphere) in liquid methane according to the 
precipitation rates of Schneider et al. (2012). 


Latitude 

min 

1ST 

(kEyrs) max (kEyrs) 

min 

RST 

(kEyrs) max (MEyrs) 

80°N 

1.6 

20.9 

3.9 

22.1 

70° N 

3.5 

45.2 

8.5 

47.8 

60°N 

7.7 

99.0 

18.7 

104.6 

60° S 

27.6 

354.9 

67.1 

375.1 

70° S 

4.1 

53.2 

10.1 

56.2 

80° S 

1.8 

23.7 

4.5 

25.0 


Table 4: Thermodynamic and physical parameters of minerals, water and ions considered 
in the dissolution reactions, collected from the literature {Lide, 2010; Langmuir, 1997; 
Ford and Williams, 2007; Brezonik and Arnold, 2011; White, 2013). 


Name 

Formula 

Mi 

[g/mol] 

Zi 

Ui 

[A] 

bi 

[kg/mol] 

AfG° 

[k J / mol] 

AiH° 

[kj / mol] 

minerals 

calcite 

CaCOs 

100.087 

0 

- 

- 

-1129.1 

-1207.6 

dolomite 

CaMg(C 03)2 

184.401 

0 

- 

- 

-2161.7 

-2324.5 

gypsum 

CaS 04 . 2 H 20 

172.171 

0 

- 

- 

-1797.2 

-2022.6 

anhydrite 

CaS 04 

136.141 

0 

- 

- 

-1321.7 

-1434.1 

halite 

NaCl 

58.443 

0 

- 

- 

-384.1 

-411.2 

aqueous species 

chlore 

Cl" 

35.460 

-1 

3.71 

0.01 

-131.3 

-167.2 

sodium 

Na+ 

22.990 

+1 

4.32 

0.06 

-261.9 

-240.1 

hydroxyl 

OH" 

17.007 

-1 

10.65 

0.21 

-157.3 

-230 

hydrogen 

H+ 

1.008 

+1 

4.78 

0.24 

0 

0 

calcium 

Ca®+ 

40.080 

+2 

4.86 

0.15 

-553.54 

-542.83 

magnesium 

Mg2+ 

24.320 

+2 

5.46 

0.22 

-454.8 

-466.8 

sulfate 

so^ 

96.060 

-2 

5 

-0.04 

-744.6 

-909.2 

carbonate 

col" 

60.010 

-2 

5.4 

0 

-527.9 

-677.1 

bicarbonate 

HC03" 

61.020 

-1 

5.4 

0 

-586.8 

-692 

carbonic acid 

H 2 CO 3 

62.025 

0 

- 

- 

-623.2 

-699.7 



solvent and 

gas 




water 

H 2 O 

18.015 

0 

- 

- 

-237.18 

-285.83 

carbon dioxide 

CO 2 

44.010 

0 

- 

- 

-394.4 

-393.51 
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Table 6: Molalities in mol/kg of solvent of minerals in water at 25° C according to the 
Ideal and the Electrolyte Solutions Theories. Pco 2 — 0-33 matm refers to normal dry air 
(Langmuir, 1997) while Pco 2 — 0-01 more refers to tropical climate environments 
(Fleurant et al, 2008), which actually could even reach considerable CO 2 content in the 
soils of up to Pco 2 = 0-11 (Ford and Williams, 2007). 


Name 

1ST 

EST 

Dissolution by dissociation 

Halite 

Gyspum 

Anhydrite 

Calcite 

Dolomite 

6.27 

8.63x10"^ 

6 . 86 x 10 "^ 

7.10x10"® 

3.77x10"® 

6.27 

2.26x10"® 

1.62x10"® 

7.72x10"® 

4.11x10"® 

Acid dissolution 

Calcite, PCO 2 = 0.33 matm 
Dolomite, PCO 2 = 0.33 matm 

4.99xl0"'‘ 

2.59x10""^ 

5.42x10-* 

2.76x10"* 

Calcite, PCO 2 = 0.01 atm 
Dolomite, PCO 2 = 0.01 atm 

1.55x10"® 

8.09x10""^ 

1.80x10"® 

8.99x10"* 

Calcite, PCO 2 = 0.11 atm 
Dolomite, PCO 2 = 0.11 atm 

3.46x10"® 

1.80x10"® 

4.28x10"® 

2.09x10"® 


Table 7: Melting and boiling points thermodynamic constants of Titanian materials. Data 
gathered from hide (2010) and Yaws (1996). The melting properties of CnHnN have been 
taken in Chirico et al. (2007). 


Formula 

Tn 

[K] 

A,^H° 
[kj / mol] 

Tb 

[Kj 

AvH° 

[k J / mol] 

CH 4 

90.680 

0.940 

111.670 

8.190 

C 2 H 6 

90.360 

0.582 

184.570 

14.690 

CsHs 

85.460 

3.500 

232.040 

19.040 

C 4 H 10 

134.850 

4.660 

272.660 

22.440 

C 2 H 2 

192.450 

4.105 

188.450 

16.674 

C 2 H 4 

104.000 

3.350 

169.380 

13.530 

C 3 H 4 

170.450 

5.349 

249.950 

22.185 

C 4 H 4 

227.600 

7.687 

278.250 

22.470 

C 4 H 6 

164.240 

7.980 

268.740 

22.470 

CeHg 

278.640 

9.870 

353.240 

30.720 

HCN 

259.860 

8.406 

298.850 

26.896 

C 2 N 2 

245.320 

8.110 

252.100 

23.330 

CH 3 CN 

229.330 

8.160 

354.800 

29.750 

C 2 H 3 CN 

189.670 

6.230 

350.450 

32.600 

C 2 H 5 CN 

207.150 

5.030 

350.500 

31.810 

C 2 H 8 N 2 

215.950 

10.070 

337.050 

28.477 

CiiHiiN 

326.650 

20.418 

538.855 

50.005 

H 2 O 

273.150 

6.010 

373.200 

40.650 

CO 2 

216.592 

9.020 

194.600 

15.326 
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Table 9: Critical properties of Titanian materials {Lide, 2010). The Rackett parame¬ 
ters (.^ra) have been gathered from the CHemical and Engineering Research Information 
Center (CHERIC) database (http://www.cheric.org/research/kdb/hcprop/cmpsrch. 
php) and from Infotherm (http://www.infotherm.com/servlet/infothermSearch). 
The compressibility factors {Zc) have been calculated for comparison with .^ra- The crit¬ 
ical properties of cyanogen and our tholins analogues (1-1 dimethyl-hydrazine C 2 H 8 N 2 , as 
suggested in Cordier et al. (2013b), based on the experiments of Quirico et al. (2008), and 
quinoline CnHuN, as sugested by Coll et al. (1995)) have been collected in Yaws (1996). 


Formula 

M 

[g/mol] 

R 

[K] 

Pc 

[MPa] 

Tc 

[cm®/mol] 

Zb.a 

^c 

CH 4 

16.043 

190.564 

4.599 

98.600 

0.289 

0.286 

C 2 H 6 

30.069 

305.320 

4.872 

145.500 

0.281 

0.279 

CsHg 

44.096 

369.830 

4.248 

203.000 

0.277 

0.280 

C 4 H 10 

58.122 

425.160 

3.787 

255.000 

0.273 

0.273 

C 2 H 2 

26.037 

308.300 

6.138 

112.200 

0.271 

0.269 

C 2 H 4 

28.053 

282.340 

5.041 

131.000 

0.282 

0.281 

C 3 H 4 

40.064 

402.000 

5.630 

163.500 

0.272 

0.275 

C 4 H 4 

52.075 

455.600 

4.860 

218.710 

0.282 

0.281 

C 4 H 6 

54.091 

425.000 

4.320 

221.000 

0.271 

0.270 

CsHs 

78.112 

562.050 

4.895 

256.000 

0.270 

0.268 

HCN 

27.026 

456.700 

5.390 

139.000 

0.197 

0.197 

C 2 N 2 

53.034 

399.900 

6.303 

144.520 

0.266 

0.274 

CH 3 CN 

41.042 

545.500 

4.850 

171.000 

0.199 

0.183 

C 2 H 3 CN 

53.063 

536.000 

4.559 

214.061 

0.261 

0.219 

C 2 H 5 CN 

55.079 

561.300 

4.260 

229.000 

0.263 

0.209 

C 2 H 8 N 2 

60.098 

522.200 

5.400 

218.241 

0.265 

0.271 

CiiHiiN 

157.212 

775.100 

3.353 

514.500 

0.247 

0.268 

H 2 O 

18.015 

647.140 

22.064 

56.000 

0.233 

0.230 

CO 2 

40.010 

304.130 

7.375 

94.000 

0.272 

0.274 


Table 10: Estimated subcooled liquid molar volumes (in cm^/mol) and associated mass 
densities for our best molar volume estimates (in g/cm^, for Eq. 24 with Zra) for Titan’s 
solids and liquids, evaluated using the Rackett equations and parameters at 91.5 K. 


Formula 

Eq. 23 
(Zra) 

Eq. 23 
(Z.) 

Eq. 24 
(Zra) 

Eq. 24 
(Zc) 

2a 

(Zra) 

P 

CH 4 

35.557 

34.928 

35.212 

34.928 

0.517 

0.456 

C 2 H 6 

46.520 

45.968 

46.229 

45.968 

0.456 

0.650 

C 3 H 8 

61.390 

62.865 

62.153 

62.865 

1.219 

0.710 

C 4 H 10 

75.882 

75.978 

75.932 

75.978 

0.080 

0.766 

C 2 H 2 

34.753 

34.185 

34.453 

34.185 

0.469 

0.756 

C 2 H 4 

42.344 

42.147 

42.240 

42.147 

0.162 

0.664 

C 3 H 4 

48.267 

49.354 

48.828 

49.354 

0.898 

0.821 

C 4 H 4 

67.280 

66.398 

66.824 

66.398 

0.729 

0.779 

C 4 H 6 

65.652 

65.175 

65.405 

65.175 

0.394 

0.827 

CgHe 

74.256 

73.269 

73.748 

73.269 

0.817 

1.059 

HCN 

30.232 

30.323 

30.279 

30.323 

0.075 

0.893 

C 2 N 2 

41.034 

43.436 

42.263 

43.436 

1.987 

1.255 

CH 3 CN 

40.219 

34.105 

36.956 

34.105 

5.041 

1.111 

C 2 H 3 CN 

71.655 

50.733 

60.016 

50.733 

17.199 

0.884 

C 2 H 5 CN 

81.023 

51.731 

64.373 

51.731 

24.046 

0.856 

C 2 H 8 N 2 

60.828 

63.522 

62.197 

63.522 

2.230 

0.966 

CiiHiiN 

123.135 

144.275 

133.476 

144.275 

17.564 

1.178 

H 2 O 

14.087 

13.692 

13.884 

13.692 

0.327 

1.298 

CO 2 

28.789 

29.225 

29.017 

29.225 

0.360 

1.379 


47 











Table 11; Subcooled liquid molar volumes (in cm^/mol) and associated mass densities (in 
g/cm^) of Titan’s solids estimated using the Rackett equation 24 with Zra at 91.5 K com¬ 
pared to the solid molar volumes (in cm^/mol) and associated mass densities (in g/cm^) 
computed from an empirical linear fit between predictions of the model and experimental 
data for some Titanian materials._ 


Name 

Formnla 

(Eq. 24, Zra) 

Rm (fit) 

2cr 



n-butane 

C 4 H 10 

75.932 

68.737 

7.195 

0.766 

0.846 

acetylene 

C 2 H 2 

34.453 

32.471 

1.982 

0.756 

0.802 

benzene 

CeHe 

73.748 

70.487 

3.261 

1.059 

1.108 

hydrogen cyanide 

HCN 

30.279 

25.562 

4.717 

0.893 

1.057 

acetonitrile 

CH 3 CN 

36.956 

37.122 

0.166 

1.111 

1.106 

water ice Ih 

H 2 O 

13.884 

19.329 

5.446 

1.298 

0.932 

carbon dioxide ice 

CO 2 

29.017 

26.109 

2.909 

1.379 

1.532 
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